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This paper discusses models of inflation based on global supersymmetry. It is shown that there 
are parameter ranges, consisent with observational constraints, for which warm inflation occurs and 
supergravity effects can be neglected. There is no need for any fine tuning of parameters. The 
thermal corrections to the inflaton potential are calculated and it is shown that they do not alter 
the warm inflationary evolution. 
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Attempts to build inflationary models based on supersymmetric Grand Unified Theories run into difficulties caused 
^ | by the size of the supergravity corrections to the inflaton potential. Inflation requires severe flatness conditions 
on the potential, but these conflict with the F-term supergravity corrections. The solutions to this problem have 
meant considering models with special cancellations or models where a different supergravity correction, the Z?-term, 
dominates (see 1] for a review). 

A totally different solution to the problem of supergravity corrections has recently been put forward, which is based 
on the realisation that the dissipation associated with warm inflation relaxes the constraints on the flatness of the 
£C) • potential |2| ■ Warm inflation can exist in a parameter regime where the supergravity corrections to the potential can 
\ be safely ignored. 

In warm inflation, particle production during inflation provides a damping effect on the inflaton. This idea has 
been around for a long time [J [3- 0- ■ but the general features of this scenario were described in H . The non- 
equlibrium dynamics has been expensively developed in subsequent work 0, and several 

^ phenomenological warm inflation models have been discussed in the literature [161 Il7l llcl ll9L |20j . 

Models in which warm inflation appears to occur spontaneousiyhave the inflaton decaying by two stages, the first 
stage into a heavy particle and the second into a light particle (2, • An example is provided by the interaction 
| Lagrangian density, 

Ph! i 

P ■ £i = ~9 2 4> 2 X 2 ~ h X^y^y (1) 

<d : 

' where x is a heavy boson and y (field tp y ) is a light fermion. The dissipation is associated with <f> — ► x ~ * VV- 
^> The ultimate destination of the vacuum energy from the inflationary phase is into excitations of the light sector 

fields. In warm inflation, one has to consider the possibility that these excitations enhance the loop corrections to 
the inflaton potential and violate the flatness conditions which inflation requires. Because of this concern, we have 
calculated the loop corrections to the potential under the assumption that the light fields thermalise. 

The light sector will typically have coupling terms representing self-interactions, or interactions with other light 
fields, in addition to the couplings given in eq Q). The relaxation time of the radiation should therefore be independent 
of the damping mechanism which is affecting the inflaton. Whether the radiation thermalises during inflation is 
therefore rather arbitrary. We assume thermalisation, but some consequences of non-thermalisation are mentioned in 
the conclusion. 

We shall consider the simplest inflationary models which include global supersymmetry. These models divide 
naturally into two classes. In the first class, which we call pure, the vacuum energy is associated only with the 
inflaton field. We find that normalising the density perturbation amplitude to the cosmic microwave background 
implies a mass scale of up to lQ ll GeV and coupling constants g and h around 0.1. 

In the second class of models, part of the vacuum energy can be linked to a false vacuum of the x held. These are 
the supersymmetric hybrid models of inflation |22l |23| . We find that normalising the density perturbation amplitude 
to the cosmic microwave background implies a mass scale of up to li) 1A GeV for the false vacuum energy and coupling 
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constants again around 0.1 . Consequently, F-tetm supersymmetric inflation with parameters in this range is of the 
warm inflationary type. 

Whilst the pressent work was nearing completion, we learned that an independent study of hybrid models of warm 
inflation was also underway, being conducted by Arjun Berera and Mar Bastero-Gil j24(. 

II. SUPERSYMMETRIC MODELS 

A. Potential and interaction terms 

We have just described how the warm inflationary scenario arises when there is a two stage reheating process 
involving a heavy boson. Global SUSY models can easily be constructed which provide the required interactions 
Consider the superpotential 

W = ~g$X 2 , (2) 

where the scalar field components of the chiral superfields and X are tp and x respectively. The scalar interaction 
terms in the theory are unpacked from the superpotential using 

c s = -\d*w\ 2 - \d x w\ 2 (3) 

We identify the inflaton with = v^Re <p, and then 

C s = -g 2 \x\ A -2g 2 4 ) 2 \x\ 2 (4) 

Supersymmetry breaking now plays an important role in determining the shape of the inflaton potential along the 
flat direction x — !]■ A 'new inflation' type of model [2^, [2(| [27j results from introducing a soft SUSY breaking 
mass M s for the \ field. The inflaton potential V{4>) is determined by the one loop correction |2*8|. 

VU>) = \?MZU^& + <&-A (5) 



2' 



Supercooled inflation requires <fro > m p , but in this parameter range the inclusion of supergravity F-term corrections 
would typically prevent the inflation from occurring. 

Hybrid inflationary models can be constructed if we change the superpotential slightly [22, H3| , 

W = 3<M 2 - g<PX 2 + g<$>X' 2 . (6) 

where A is a constant and X and X' are a pair of superfields. The interaction terms are now 

Cs = -g 2 \ X 2 ~ X' 2 ~ AT - 2g 2 tf(\x\ 2 + \X'\ 2 ) (7) 

In hybrid models, the x field is stable at % = during inflation and the potential is dominated by the constant term 
<7 2 A 4 . The x field becomes unstable at when <fi falls below the critical value 4> c = A. 

The supersymmetry is broken by the inlaton field resulting in a non- vanishing one loop contribution to the inflaton 
potential. The presence of the second superfield helps produce a potential which is suitable for inflation by reducing 
the size of the quantum corrections. For <j) 3> A, 

y^) = 9 2 ^ + ^4 2 4f] (8) 



4tt 2 V A 2 

The heavy sector plays a double role in contributing to the vacuum energy and damping the inflaton field. |44j| 

For an efficient two stage reheating process, we introduce an additional light sector Y, which can be coupled through 
a superpotential 

W = -g$X 2 - hXY 2 (9) 
The Yukawa interaction terms are recovered from 

1 d 2 W — 1 d 2 W* — 

2 d(p n d(j) m 2 d4>* n d(j)* m 

where <f> m is a superfield and Pi, = 1 — Pr = (1 +7g)/2. The interactions contain terms such as those in eqn. (JIJ, and 
lead to a friction term oc (f> in the inflaton field equation . They also have an effect on the vacuum polarization of the 
X field, which in turn can affect the inflaton potential. The full set of interaction terms and the vacuum polarisation 
are discussed in section 3. 
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B. Inflationary dynamics 



In an expanding, homogeneous universe, the inflaton equation of motion is given by 

4>+(3H + T)j> + V Tt<j> =0, (11) 

where Vr(<p, T) is the thermodynamic potential and T(<fi,T) is the damping term due to interactions between the 
inflaton <j> and surrounding fields. For supercooled inflation, this damping term is negligible compared to the Hubble 
damping term. The interesting regime of warm inflation is characterised by large damping terms. To distinguish 
between the two inflationary scenarios, a dimensionless parameter, r, is introduced to denote the relative strength 
between the damping terms 

r 

r ~ W 

We shall take warm inflation in the limit r>l. 

As the inflaton evolves, energy dissipates into radiation and entropy is produced. Simple thermodynamic relations 
lead to a definition of entropy density, 

4tt 2 

s(d>,T) = -V T , T = — 5 *T 3 + ..., (12) 

where <?* is the effective particle number and the dots denote contributions from the thermal correction to the potential. 
In the warm inflationary scenario, inflation is characterised by three slow-roll equations |2^| 

1 _ V T,4> rn„ _ !2 



3iJ(l 



Ts = r<j)\ 3H z = 8nGV T . (13) 



The second equation denotes conservation of energy, while the third is the usual Friedmann equation. Slow-roll 
automatically implies inflation, a > 0, and the consistency of slow-roll is governed by a set of slow-roll parameters: 



1 ' V tA m _ml(V T ^\ ? _ml / r,^ \ 5 _ TVtjt ^ (m) 



€ ~ 16tt V V T J ' '' 8tt V V T J ' " 8tt V TV T J ' " V T . 

where m~ 2 is Newton's constant. The slow-roll approximation is consistent when the above par ameters are less than 
r. Supergravity F-term corrections, without special cancelations, lead to r\ of order unity |23. We shall therefore 
concentrate on the range, 

1 < r\ < r. (15) 

The thermal corrections to the potential will be calculated in the next section. For the remainder of this section we 
shall examine the situation where the corrections are small, Vr — V and the slow-roll parameter 6 = 0. 

An important observational constraint on the model is set by the density perturbation amplitude. In our case, where 
we have assumed that the radiation has thermalised, the thermal fluctuations induce scalar density fluctuations. The 
amplitude A can be obtained analytically [2||[3(j, and for r ^> 1, 

i/l/ 4 -3/4 1/4*2/3 fins 

V = ar h e h A m p ( 16 ) 

1/12 

where a w 0.68g* and the parameters are evaluated at the time th that the perturbation scale crossed the horizon. 
The value of A infered from cosmic microwave observations is around A i=s 5.4 x 10~ 5 on the 500 MPc scale [3l|. 

Limits on the mass parameters can be found by combining the slow-roll limits l|15l) with the constraint from the 
density perturbations (|16fl . The pure inflation model with potential (JSJ has two mass parameters g M s and </>o- Order 
of magnitude estimates can be obtained by taking 0^ ~ (f> for the value of 4> at horizon crossing (which is consistent 
with numerical solutions |28||1 and e w r\. The normalisation condition l|16|) gives 

gM s » 3.3 x 10- 6 rj h r- 3/2 m p . (17) 

The upper limit for gM s set by ljT5)l is of the order 10 13 GeV. 

For hybrid inflation, a similar approximation can be made when the vacuum energy g 2 A 4 dominates the potential 
and 4> ^> A. In this case, 

5 2 ff 4 m 2 
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The normalisation condition l|16(l gives 



A m 3.88 x varSi!V ,i m p ( 19 ) 



When combined with the conditions > T)h > 1, the upper limit on A is of order 10 15 GeV. 

More detailed limits can be placed on the parameters when we know the form of the friction term in the inflaton 
equation. For the interactions in eqn. l(TJl. the friction term has been calculated in the zero-temperature limit |2l) 
and is given to leading order in h by 

r « 70- (20) 
The value of 7 depends on the decay process. For % — > 2y and x 2y, the leading order contributions to 7 are 

7 (x-2y)=7(x-2y) = ^|!^ (21) 



For the fermionic channel \ VVi 



Hence the total 



7 16tt 2 

To be in the perturbative regime, with g < 1 and h < 1, sets a requirement 7 < 8.9 x 10~ 3 . 
We can relate r to 7 using the slow- roll equations (|14|l . 

r 7^ to p 

r = — - — 



(23) 



3H (24TTV) 1 / 2 



(24) 



The non-hybrid models can be regarded as having three parameters gM s , 0o and 7- The normalisation from the 
density fluctuation amplitude provides one constraint which can be used to eliminate one parameter, let's say </>o- The 
warm inflationary regime can then be displayed as bounds on the remaining two parameters. 

The horizon crossing timescale depends on the number N of e-folds of the scale factor before the end of inflation. 
As a rough guide, we can take = Nrjh- The normalisation condition (|16|l gives 

0o = 8.2 x 10 _6 7 _1 N-V 2 m p (25) 

Note that 0o < m-p is needed for rj to remain larger than one throughout the inflationary era, which is a necessary 
requirement for the neglect of supergravity corrections. After eliminating </>o, there is a consistency requirement 

7 gM a < 4.8 x lQ- ll N- 2 m p (26) 

for warm inflation. 

The hybrid models can be taken to have parameters M — g x ^ 2 K, g and 7. The horizon crossing timescale is given 
approximately by the relation r% = ZNrjh- The normalisation condition Hl(jfl allows us to express g as 

g w 1.5 x 10 21 7 2 A/ 2 7V 5/2 m~ 2 (27) 

The condition 77 > 1 becomes 

7 > 1.3 x 10- e N~ 1/2 (28) 

Taken together with the expression for 7 given by eq. (|23|l . this gives lower limits on g and h. For example, if g « h, 
then g > 0.1. Another condition follows from (f> > A, which imples a consistency condition 

7M > 3.8 x 10" 12 iV- 13/8 mp (29) 

for the warm inflationary regime. 
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FIG. 1: The figures show, approximately, the allowed values (region A) of the parameters 7, gM a and g 1 ^ 2 A for (a) pure and 
(b) hybrid inflation. The line SUGR shows the limit for which SUGR corrections can be neglected. In case (a), the limit of 
consistent warm inflation is shown as WI. In case (b), the limits for consistent warm inflation are <f>h = A (WI1) and g = 1 
(WI2). 



The parameter ranges are shown in figure ^ The number of e- folds of inflation has been taken to be TV = 60. It is 
clear from the figure, that the warm inflation occurs for a broad range of parameters. In the case of hybrid inflation, 
range includes 0.03 < g < 1 which means that many F-tevm inflationary models should be have warm rather than 
cold inflation |23j. 

Finally, we can use the slow-roll equations to determine what the temperature is during the inflationary phase. The 
slow-roll equations (fH5|) give 

Ts = \r- x eV (30) 
Once again making use of the normalisation provided by the density fluctuation amplitude, and the slow-roll relation 

v - 8wTlV rm 

V >** ~ ~mf> (31) 

we can compare Th to the slope of the potential, 

T h ^230r 1 h /2 r 1 - 1/2 (V 4>4 ,) 1 ^ (32) 

Note that Th is always larger than the mass scale responsible for the slope of the potential. In this situtation, wc 
should be concerned that thermal corrections to the potential may make the models untenable. The calculation of 
these thermal corrections to the potential is therefore necessary. 

In the case of pure inflation, V.^ ~ 2g 2 M 2 and the temperature during inflation is approximately 

T h ~ 320 N 1/2 gM s . (33) 
In the case of hybrid inflation, using eq. (|30J) and eq. Ijl8|l , 

T h ^f^y\^N-^M (34) 



Note that the vacuum energy of the hybrid model after inflation, before the second field decays, is approximately 
M 4 . The values of the temperature are sufficiently high such that, if we have local s upersym metry, then we are in 
danger of violating constraints set by the thermal producation of gravitinos |32l I3H , I34I l35l l36l l37j| . These constraints 
can be satisfied by taking small values of the mass parameters, corresponding to large values of r 38] . The allowed 
parameter ranges permit this, but so far this appears to be an unnatural feature of the models under discusion. 
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III. THERMAL CORRECTIONS 

The full set of interaction terms obtained from the superpotential © are 

Cs = ,9 2 (A 2 - |x| 2 ) 2 + VM 2 M 2 + 4^y| 2 | X | 2 + h 2 \yf + 2 5 %Vx t + V^VX) (35) 

+2 9 (x^ x PltP v + xH x Pr1> v ) + 2h(y%P L iJj x + y^ y P R ^ x ). (36) 

Thermalisation conditions of the light fermion, y depend directly on the mass of the fermion and its self-interaction. 
These properties are inherent to quadratic and cubic terms in the superpotential (n x Y 2 , Ay 3 ), which we have 
not specified. We assume that the interactions are such that the light fermions thcrmalise and we calculate the 
corresponding thermal effects. The ip y interactions to fields other than the inflaton, <fi or x fields, will have no effect 
on the thermal corrections to the inflaton effective potential. We can therefore disregard the exact nature of these 
interactions. 

For the model considered, if y thermalises but x an d V'x ^° n °tj then important simplifications can be made. 
The thermal corrections in the action appear as a result of the self-energies of the X an d X fields. Inside the self- 
energy loops, y and y are taken to be very light, so the Hard Thermal Loop (HTL) approximation can be made, i.e. 
T » m y , rtiy. Outside the loop, however, the x fields are heavy, and T <C m x , m x . 

We use the imaginary time formalism and adopt the notation that 4-momenta are written in upper case and 3- 
momenta as written in lower case, so that P M = (u), p). The boson and fermion propagators of the x fields are G and 
S respectively. The contribution to the effective potential of the inflaton field from the x fields is given by 

V * = lwf ln dct ( G ") -/ J^f ln clet^S*- 1 ) 172 (37) 

after regularisation has been applied. The detailed calculation of the thermal corrections to the fermionic and bosonic 
masses m x and m x will be given in the following sections. 

A. Fermion Contribution 




i) x <ip y Vx 



FIG. 2: Diagram for the fermionic self energy 



If we set y = (yi + iy-2)/y2, then the yi"4> y ip x terms in eqn. (|36[) lead to two similar fermion self-energy diagrams 

(fig. EJ) for x- The vertex factors are iy/2h and \ 2hrf§ respectively. Using the properties of 75, contributions from both 
diagrams are found to be identical. Thermal feynman rules applied to the diagram in figure [21 lead to an expression 
for the fermion self-energy, E, 

{ ri 3 k ~ 
E(P) = -^ T ll J J^W-P)^{K)^{P-K). (38) 

where A(if) « K~ 2 , k° = 2nnT for bosons and k° = (2n + 1)ttT for fermions (denoted by a tilde). 

Many similar diagrams appear in the literature, for example [39l liol EllE^ . specifically with regards to the HTL 
loop corrections to electron propagators in QCD. The calculation given here follows |4l| with only slight changes, due 
to having a scalar field rather than a vector field. Therefore we only need quote the result, 

2p \p J 2p 







1 Qo 


\~p) 


P 





(39) 
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Note that the overall factor of 1/2, which is different to the literature, is a convention we have adopted for the 
definition of rrij for later convenience. Qq(x) is the Legendre function of the second kind, 



Qo{x) 



1, ar + 1 
- In 

2 at-1 



Accounting for the contribution of two diagrams, the fermion thermal mass is 

, h 2 T 2 



In our conventions, the inverse propagator is iS 1 = f— m% — £ and thus can be written 



where 



iS 1 = A 7o - A s j ■ p- nix 



1 f ^ I lUJ 
A = iuj - —Qo 
2p 



2p 



Hence the combination 



(SS*)- 1 = -Al + A 2 s +m\ 



The fermionic contribution to the effective potential (|37[1 becomes 

d 4 P 



V f = -2 



(27T)' 



In 



my im 



We can obtain the leading terms using m x 3> rrif , 

f d 4 P 



{2*y 



\n{oj 2 +p 2 +m| +m)} + 0(mj) 



The evaluation of this regularised integral is a standard excercise 

f 



/ 



32tt 



2 (m\ + 2m^m 2 ) In ( -* 



f 



where fi has been introduced by the regularisation. 



(40) 

(41) 

(42) 

(43) 
(44) 

(45) 
(46) 

(47) 
(48) 



B. Boson Contribution 



Three diagrams are expected to contribute to the bosonic self-energy as shown in figure [3| Due to the interaction 
terms Xi^Py^Py m ec L n - IpSt ; there will be two diagrams similar to figure 3(a) with vertex factors iV2h and \/2hrf§ 
respectively. The two diagrams result in identical expre ssions, each with a symmetry factor of 1/2 because the fermions 
are Majorana. The self-energy for xi from figure 3(a) is given by 



In the HTL limit, 



H(P) a ^h'T^J j^Tr \m- f)\ A(K)A(K P) 
n(P) a = -Ah 2 TJ2J j^K 2 A(K)A(K P) 



(49) 



(50) 
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(b) 

FIG. 3: Diagrams for the bosonic self energy 
Using A(K) ~ K~ 2 the final result is a single fermion loop expression, 

3i 




(51) 



The second contribution stems from the four tadpole diagrams coming from the terms XiUi m ecm - 1|35|) . Each 
diagram is identical, and has a symmetry factor of 1/2. The self energy is given by the expression, 

d 3 k 



(52) 



The Feynman diagram in figure 3(c) does not result in an 0(T ) contribution to the self-energy, but instead gives 



a possibly important contribution of 0(m 2 ). The self-energy contribution is given by 

U(P) C = Ag 2 h 2 4> 2 T^2j -0pA(K)A(K P) (53) 



(27T) 

For the evaluation of this integral in the HTL limit see, for example, The result is 

„2^2 f dk 1 f k-p \ 1 

(2tt) 2 p \k + p)e kT -l 

The integral evaluuates to 



n(P) c = V/jV 2 / z ln ( irrz ) zw^r- ( 54 ) 



n(P) c » ^a^fcVlog^- (55) 

in the HTL limit. 

The inverse boson propagator is defined as G _1 = P 2 + m 2 , + II. The total contribution from the first two diagrams 
in figure |3 defines a contribution m 2 to the mass, 

m 2 = U a + U b = h 2 T 2 (56) 

Note that m 2 = 4<? 2 2 , and the temperature dependent part of II C can be regarded as a contribution to the finite 
temperature coupling constant g(T), 

.g 2 (T)=.g 2 (l + ^ 2 logJ). (57) 

Since the inflationary dynamics is not sensitive to the precise value of g, we shall not distinguish between g and g(T). 
The bosonic contribution to the effective potential is obtained from the two fields Xii 



J ^ ln { p2 + m x+™'} 



32tt 



(2tt) 

— K + m 2 ) In (58) 
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Effective Potential 



The total effective potential is given by the sum of both the fermionic and bosonic contributions. The largest of 
the 0(T 2 ) terms cancel due to the fact that mj, = rrif. This is, of course, due to the underlying supersymmetry, but 
it happens despite the fact that supersymmetry is broken at non-zero temperatures. The remaining 0(T 2 ) terms are 
due to other sources of SUSY breaking. For the case of soft SUSY breaking, we can consider the \ boson and \ 
fcrmion masses, 



TO 1 



■bf,/ + M s 
2g 2 4?. 



The leading order terms in the potential are 

1 [, o 



32tt 2 



to?) 2 In 



to 2 + to 2 



i 2 f) 2 In 



to 2 +m 2 



constant. 



(59) 
(60) 



(61) 



For our scenario, g<p 3> M' s , which results in 

V x « \m 2 s (g 2 <p 2 + \h 2 T 2 



In 



2^,2 



9 <t> 



\h 2 T* 



1 ) + -g 2 M 2 s <&. 



(62) 



We have defined Af 2 = M^ 2 /(87r 2 ), and 4>o has been chosen such that both the potential and its derivative vanish at 
<p = (f>o when T — 0. 

The thermodynamic potential of the inflaton in the models under consideration is determined primarily by the x 
loop contribution calculated above and the free energy of the light radiation fields, 



U(0,T) = -- 3 *T 4 + U x (0,T) 



(63) 



Figure 0| graphically shows the temperature dependance of the potential. The scale of the temperature corrections is 
set by Tq, where hTo = 2g<j)Q. 




FIG. 4: The thermodynamic potential of the pure inflationary model is depicted. On the vertical axis, V = Vr((f>,T) — Vr(0, T) 
and Vo = Vt(0,0). The critical temperature hTo = 2g4>o- 
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The Hybrid models © can be treated in a similar manner if we couple the two heavy superfields to the light 
superfield with equal couplings h (which is consistent with an underlying U{\) symmetry). The masses of the x an d 
X 1 fields are 

m 2 , = 2g 2 (j) 2 + 2g 2 K 2 (64) 
m 2 ,, = 2g 2 <j) 2 -2g 2 K 2 (65) 
TO | = m |, = 2g 2 (j) 2 (66) 



1 4 . 4] / 2gV + \h^ \ 
V x = 7-2 9 A In — (67) 



The one loop correction to the inflaton potential is now 

— ^-o 4 A 4 lu . 
4tt 2 V 5 A 

for > A. 

D. Friction term 

The zero temperature friction term calculated by Berera and Ramos |'2ll can be modified for our interaction terms 
with the decays cf> —* \ — > 2y and (f> — > x — > 2y, 

r : z ; ^ ^ (68) 



2tt K + T2) 1/2 (2m x (m| + r 2 ) 1/2 + 2m| 



where T x is the x decay width, 



r v=T7-"'x ( '-— fl +T77Z'"x ( 69 ) 




2^,2 



This reduces to the expression used earlier when m y = and m 2 , = 2g 

One would expect that the thermal corrections to T will manifest themselves as corrections to m x and T x . Since 
the latter already contains factors of h, the thermal corrections to T will be of order h 3 T 2 . The correction to T due 
to m x will be of order h 4 T 2 . The corrections to the effective potential are order h 2 T 2 and are therefore taken to be 
the dominant effect. Corrections to the friction term can therefore be ignored. 



E. Effects on the inflationary dynamics 

The temperature range relevant for the warm inflationary scenario was given in eq 132f) . In the models under 
discussion, h 2 T 2 <§C g 2 cj) 2 and the thermal corrections make only a small change to the height of the thermodynamic 
potential Vt- The effect on the slope of the potential is more delicate, however, and has to be investigated seperately. 
The slope can be quantified by the slow roll parameter r\. Consider the change in 77, 

V V,<1>4> 

In the non-hybrid case, 

(7D 

Vk g n 

Given eqs. (|33|l . (|25|l and l|26|) . this correction is of order 10~ 5 h 2 g~ 2 at most. 
For hybrid inflation, using l|67() . 



8y h „ h 2 T 2 
Vh g 2 4> 2 h 



(72) 



Proceeding as above, using eq. (|34f> . this correction is of order 10 7 h 2 g 2 
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IV. CONCLUSION 



We have found that the inflation occurs naturally in particle models with global supersymmetry when the dissipative 
effects of particle production are taken into account. The warm inflationary scenario escapes the flatness problems 
which arise when supercooled inflation is combined with global supersymmetry. The parameter restrictions on the 
model are not severe, with the possible exception of a gravitino constraint, and there is a correspondence between 
mass parameters required for the observed density fluctuation amplitude and the parameter values of interest for 
supersymmetric Grand Unified Theories. 

We have demonstrated that, in a two stage reheating process, the thermal corrections to the inflaton potential are 
small, due to fermion-boson cancellations . The assumptions used for the models have been relatively mild, consisting 
mainly of the following: 

• At least one superfield has vanishing, or very weak coupling, to the inflaton and another has non-vanishing 
coupling. During inflation, the former will naturally become a 'light' sector, and the latter a 'heavy' sector. 

• There is either (a) soft SUSY breaking in the heavy sector, which we called the pure inflation model or (b) a 
false vacuum energy, and two equally coupled heavy superfields, which we identify as the hybrid model. 

• We have assumed that the light radiation thermalises. 

The last assumption was needed to avoid far from equilibrium calculations. However, in the absence of thermalisation, 
the energy will still be dumped into the light sector. The important features of the light particle distribution can 
be described in terms of non-thermal occupation numbers n(k). If the boson and fermion occupation numbers are 
similar, we might still expect the cancellation of correction to the effective potential which we have found here. 

We have made use of the density fluctuation amplitude when setting limits on the parameters in the models. If 
the radiation does not thermalise, we would expect to find changes in the predicted value of the density fluctuation 
amplitude and this remains to be investigated further. 



We are grateful to Arjun Berera for discussions about this work and to Mar Bastero-Gil for pointing out ref. [23| 
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